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The use of quantum entanglement to study condensed matter systems has been flourishing in 
critical systems and topological phases. Additionally, using real-space entanglement entropies and 
entanglement spectra one can characterize localized and delocalized phases of disordered fermion 
systems. Here we instead propose the momentum-space entanglement spectrum as a means of 
characterizing disordered models. We show that localization in ID arises from the momentum space 
entanglement between left and right movers and illustrate our methods using explicit models with 
spatially correlated disorder that exhibit phases which avoid complete Anderson localization. The 
momentum space entanglement spectrum clearly reveals the location of delocalized states in the 
energy spectrum and can be used as a signature of the phase transition between a delocalized and 
localized phase. 



New light has been shed on quantum many-body 
ground states via their entanglement properties. 
Entanglement measures, such as the topological 
entanglement entropy [THE] and the entanglement 
spectrum [6HT2], with their origin in quantum in- 
formation theory have begun to uncover subtle fea- 
tures of topological insulators and systems with 
topological order. One field where entanglement 
characteristics have begun to make progress is in 
disordered materials. For example, Ref. [T3land[T4l 
showed a connection between the position-space 
entanglement entropy and a multi-fractal scaling 
exponent which was extended by Ref. [15] to estab- 
lish a connection between Renyi entropies and the 
universal part of the multi-fractal spectrum. Ad- 
ditionally Ref. [16] and [17] showed that the level- 
spacing statistics of the entanglement spectrum 
match the same statistical ensemble as the energy 
spectrum. Since entanglement in position-space 
captures information about long-range correlations 
it is natural that it should reveal information about 
localized and de-localized states. 

In this article we also consider the entanglement 
characterization of disordered fermion models, but 
from the perspective of momentum space instead of 
position space. The advantage (and from another 
perspective the disadvantage) of entanglement is 
that it crucially depends on the way Hilbert space 
is partitioned into two sectors. The partitioning 
can be done in any way, though typically the par- 
tition (or cut for short) is made in position-space. 
While the dependence on a cut makes entangle- 
ment non-unique, the benefit is that there is a 
possibility to learn more information about a sys- 
tem by taking different cuts. The notion of cut- 
ting the system in momentum space was studied 
for translationally invariant ID spin-chains in Ref. 
[SJ but, as we will mention again below, this type of 
cut yields no information at all for translationally- 
invariant, non-interacting fermion systems. While 



a momentum-space cut is useless when translation 
symmetry is preserved, we show, however, that it 
does lead to valuable information if the system is 
disordered. The disorder couples/entangles differ- 
ent momentum states through scattering and the 
momentum space entanglement spectrum shows 
clean signatures to characterize regions of local- 
ized and de- localized states. In this work we focus 
solely on ID models where the essential idea is that 
right-moving states, which were delocalized in the 
translationally invariant limit, remain delocalized 
if the disorder does not entangle them with left- 
movers. 

To illustrate localization transitions in ID we 
must use special classes of models as it is well 
known that disordered, ID electronic systems be- 
come insulating due to the phenomenon of Ander- 
son localization. An exception to Anderson lo- 
calization in ID was discovered by Dunlap et al. 
who predicted that a model with spatially corre- 
lated disorder, the so-called random dimer model 
(RDM), could avoid complete localization (at least 
in finite-size systems) [18] . The RDM model is de- 
scribed by a tight-binding lattice Hamiltonian 

N N 

H = Yl 1 ( C Ul C n + C n C n+l) + Yl 6 n4i C n, (1) 
n=l n=l 

where n labels the N sites of the lattice, t is the 
nearest neighbor hopping , and e n represents a site- 
dependent energy. In the RDM the site energies, 
which represent the random disorder have a spe- 
cific structure: (i) e n can take two possible values 
e a or 65 (ii) the two energies are randomly placed 
on each site with the constraint that when e a is 
chosen it always placed on two consecutive sites 
(thus a dimer, see Fig. [3J. This constraint is what 
renders the disorder spatially correlated. Dunlap 
et al. showed that the RDM model possesses a 
pair of delocalized eigenstates that can propagate 
throughout the lattice without being localized by 
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FIG. 1. Entanglement spectrum as the Fermi energy 
is varied, with (a)e& = O.Ot, (b) — O.lt, and (c)— 3.0£. 
The red dashed lines denote the resonant energy, (d) 
Entanglement spectrum as the disorder strength is var- 
ied, with Ef = fixed at the resonant energy. The red 
dashed lines denote the critical point beyond which the 
extended states of the random dimer model are lost. 



the disorder. These states are degenerate and oc- 
cur at an energy E = e a , which is termed the res- 
onant energy. An additional 0(y/~N) states in the 
neighborhood of the resonant energy remain ex- 
tended as well, and as a result, a finite fraction 
of the eigenstates are able to resist Anderson lo- 
calization. Thus, a measurable conductance can 
in principle be observed in finite-length wires or 
polymers [19]. 

The RDM can be further generalized to the "ri- 
mer" case in which chains of length n and en- 
ergy e a are randomly placed throughout the lat- 
tice; the n = 1 case corresponds to uncorrelated 
disorder, and n = 2 to the RDM. The general 
case has n — 1 resonant energies at E n (m) = 
e a — 2t cos (irm/n) , m = 1, • • • , n — 1. In what 
follows, we will mainly focus on the n — 2 dimer 
and n = 3 trimer (RTM) cases. Without loss of 
generality we set e a = so that the resonance en- 
ergy £^2 = (E3 = d=t) for the dimer (trimer) 
model and the disorder strength is quantified by 
£b — £a — £b- For the RDM the delocalized states ex- 
ist for 1 65 1 < 2t and for the RTM delocalized states 
at E 3 = -t(+t) for -3t < e b < t(-t < e b < 3t). 

In order to make contact with previous work, 
and for comparison to our main results, we will 
first consider the spatial entanglement of the RDM. 
We divide the Hilbert space into two subspaces A 
and S, which contain sites confined to the region 
n e [I, N/2] and n e [N/2 + l,iV], respectively. 



We will focus on the single-particle entanglement 
spectrum calculated via Peschel's method which 
considers the eigenvalues Q of the correlation ma- 
trix Cij = (Q\c[cj\Q) where \Q) is a free-fermion 
ground state and i,j e [1, N/2] [20 . The {d} form 
the single-particle entanglement spectrum and al- 
ways lie between and 1 with eigenvalues near 1/2 
yielding the most entanglement. When we calcu- 
late momentum- space entanglement later the only 
difference is that z, j in dj will refer to different 
single-particle momenta, which are discrete indices 
because TV is finite, and they will be restricted to 
lie in certain regions of momentum space instead 
of position space. 

In Figs. [TJi,b,c we show the entanglement spec- 
trum as a function of the Fermi energy for disor- 
der strengths fixed at the values e b = 0.0£, —O.lt, 
and — 3.0t respectively (note that only one disor- 
der realization is shown and no averaging has been 
performed). For the clean system we see large en- 
tanglement for all values of the Fermi-level within 
the band. For e b = —O.lt we still see a lot of en- 
tanglement and one can see by eye that there is 
level-repulsion in the entanglement spectrum for 
Fermi-energies in a window around the resonant 
energy E = 0. This is reminiscent of the level re- 
pulsion commonly seen in the energy spectrum for 
energy ranges containing delocalized states, and 
which was confirmed to exist in the entanglement 
spectrum in Ref. [16 . Finally, for e b = — 3.0t, 
which is tuned to the localized insulating phase, we 
see that the entanglement is suppressed compared 
to the other figures with eigenvalues clustered near 
C = and C = 1 which heuristically signals the ab- 
sence of delocalized states. In Fig. [T]i, we show the 
entanglement spectrum as a function of disorder, 
with the Fermi level tuned exactly to the resonance 
energy. Spatial entanglement appears to be more 
pronounced when |e&| < 2t, as expected in the de- 
localized phase. There is a somewhat clear distinc- 
tion between the localized and delocalized phases 
based on the density of entanglement eigenvalues 
away from £ = 0, 1. In fact, even the highly entan- 
gled modes in the localized phase are most likely 
due to local entanglement i.e. localized states near 
the cut which happen to have weight on both sides 
of the partition. 

While some features of the phases of the RDM 
are apparent from the spatial entanglement, we 
will show that performing a momentum-space cut 
leads to more explicit results. A trivial, but ulti- 
mately illustrative example is to consider a mo- 
mentum space cut in the translationally invari- 
ant limit. In this limit, since momentum is con- 
served, entanglement is completely absent and the 
set of {Ci} will have a number of 1 (0) eigenval- 
ues equal to the number of occupied (unoccupied) 
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FIG. 2. Entanglement spectrum as the Fermi en- 
ergy is varied for (a) disordered dimers (b) disordered 
trimers for eb — — 0.5£. Entanglement spectrum as the 
disorder strength is varied, for (c) disordered dimers 
(d) and disordered trimers for 22/ at resonance. 



single-particle momentum states which depends on 
the Femi level. Thus, for a clean system, where 
all states can be considered as delocalized plane- 
waves, there is no entanglement. Once disorder is 
turned on, momentum states are mixed and the 
{Q} will have more generic eigenvalues, but the 
key to locating regions of extended states is to 
search for Fermi-levels where the entanglement is 
suppressed. This is the opposite of the spatial en- 
tanglement paradigm. The other important con- 
cept is that the clearest signatures are seen if the 
momentum cut is chosen so that region A consists 
of all the right-movers and B contains all the left- 
movers. An essential ingredient for determining 
the cut is thus the translationally invariant band- 
structure which, as we will see in detail later, is cru- 
cial for determining the the single-particle group 
velocities and thus the best momentum cut. 

For the RDM and RTM the relevant partitioned 
subspaces have momenta restricted so that kj, G 
[0, 7r] and k& G [7r, 2tt] since these regions make 
up the left and right movers respectively. In Figs. 
|2^i,b we show the resulting entanglement spectrum 
for the RDM and RTM respectively as a function 
of the Fermi level, for this choice of partition. The 
disorder strength was set at 65 = — 0M for both the 
RDM and RTM. Remarkably, the resonant ener- 
gies are revealed with great clarity where a marked 
suppression of entanglement occurs. In the vicinity 
of the resonant energies 0(y/~N) of the Q approach 
either or 1. On the other hand, for Fermi ener- 
gies away from the resonances, the spectrum fills 



up with entangled eigenvalues due to scattering be- 
tween left- movers and right- movers. In Figs. (2j},d 
we show the entanglement spectrum as a function 
of disorder strength for the RDM at Fermi-level 
E F = E 2 (l) = and the RTM at Fermi-level 
E F — 223(1) = —t respectively. It is clear from 
Fig. ^ that the entanglement is strongly sup- 
pressed in the RDM up to the phase transition 
at 1 65 1 = 2t, where the extended states are lost. 
Similarly, for the RTM we see in Fig. [2]i that the 
state at 223(1) = —t only remains delocalized for 
—3t < 65 <t. One other feature of note is the sup- 
pressed entanglement in the RTM for 65 < St. 
This effect is not generic and is due to a pathology 
of n-mer models (for n > 2) in the strong disorder 
regime, which is dominated by bound states local- 
ized on the n-mers themselves near the resonant 
energies. The position of the resonant Fermi- level 
with respect to the energy of the n discrete modes 
on each n-mer leads to this spurious suppression 
which we will not discuss any further. 

To understand this behavior of the RDM we 
study the Hamiltonian in the momentum- space ba- 
sis replacing c n = ^^Z k e zkn c k in Eq.|Tj). This 
leads to the expression 

H = ^2 £(k)c{c k + ^2 ^Afe4 c fe' ' ( 2 ) 

k k, k' 

where £(k) = e a + 2tcosk. The matrix 
elements CAfe of the disorder potential read 
e~Afc = (o> - e ) [6k,k f ~ f (Afc, {n}) S(Ak)} , where 
f(M,{n}) = -^Y,iZ iri{k - k '\ and S(Ak) = 
^= (l + e^*-*')) . Here, Ak = k f - jfe, the n 
denote the random positions of the dimers, and 
we have restored e a to show the dependence ex- 
plicitly. The function /(Afc, {r^}) captures the 
random part of the disorder since it depends ex- 
plicitly on the positions of the dimers whereas 
the structure-factor S(Ak) is a coherent sum of 
phases, independent of any randomness. The func- 
tion S(Ak) ~ (Ak — Q) + • • • , when expanded 
around Q = 7r, which, even after convolution 
with /(Afc, {ri}), suppressed the scattering events 
k — )► k' + Q as shown in Fig. [3) The gener- 
alization to the arbitrary n-mer case can readily 
be made by replacing the structure factor with 

s n (Ak) = ^Er = o etmM) which has zer ° es 

at Qn(rn) = 2nm/n, m = 1, • • • , n — 1. In Fig. [3j 
we show examples of the disorder distribution in 
both position and momentum space respectively, 
for the n = 1, 2, 3 cases. 

For generic weak disorder, scattering predomi- 
nantly mixes degenerate single-particle momentum 
states 1(71,2) • If the ^1,2) have opposite group ve- 
locity then this mixing leads to backscattering, lo- 
calization, and entanglement between left-movers 
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FIG. 3. Profile of position (left) and Fourier (right) 
distributions of the n-mer disorder potential, with n = 
1, 2, 3 (top to bottom). 



and right-movers. If we consider the RDM model, 
then states satisfying q\ — q^ = Q are never mixed. 
A resultant suppression of momentum entangle- 
ment occurs when this condition is met by degen- 
erate states with opposite group velocity. For the 
RDM the momenta at which the condition £(qi) = 
£(Q2) = £(Qi — Q) is satisfied are given by q\^ = 
±7r/2. At these momenta the resonant energy for 
the RDM is E 2 (l) =£ (qi) = e a + 2tcos (tt/2) = e a 
which is exactly what we noted earlier, though we 
chose e a = for convenience. It then follows that 
the suppression of left and right-mover entangle- 
ment matches the existence of extended states in 
the RDM model near the resonant energy. This 
connection constitutes the main result of this work. 
By analogy with the dimer case, the resonant en- 
ergies of the n-mer model are determined by sat- 
isfying £(q) = S(q — Q n (m)). Thus, we see that, 
whenever Ep is at one of the resonant energies, 
the absence of left and right-moving entanglement 
is correlated with the presence of extended states 
in the system, just like in the particular n = 2 
dimer case. 

Finally, we consider an artificial model in order 
to verify that: (i) the presence of zeroes in the 
momentum-space matrix elements generically lead 
to resonances, regardless of the model-dependent 
shape that the disorder might have; (ii) the in- 
terpretation in terms of left and right-mover en- 
tanglement and momentum space partitioning is 
consistent. To accomplish (i), we construct a new 
disorder distribution V n such that Va& = for all 
Afc, except those in some finite ranges for which we 
set the magnitude to be a nonzero, constant value. 
Disorder is introduced by randomizing the phases 
of the scattering elements in momentum space. An 
example profile of such a V n and |Va/c| are shown 
in Figs. [4^l, b for a particular realization of the dis- 
order. To address (ii), we implement an exclusively 
next-nearest neighbor hopping model, so that the 
group velocity becomes v g (k) = — 4tsin2/c. In do- 
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FIG. 4. (a) and (b): spatial and Fourier components of 
V respectively, (c) and (d): entanglement spectrum as 
a function of the Fermi level. The dispersion relation 
is overlaid on the entanglement spectrum as a guide, 
with the restrictions on A schematically emphasized 
by the solid orange segments. The dashed green seg- 
ments correspond to region £>, which is traced out. The 
black dashed lines delineate the regions over which one 
expects suppression of entanglement according to the 
zeroes of Vau- 



ing this, we consider two possible partitions of the 
Hilbert space: one choice corresponds to the re- 
striction kj± G [0,7r], which is what we chose for 
the RDM, and the other takes into account the new 
location (s) of the left movers in momentum space, 
so that kj, G [0, tt/2) U [it, Stt/2]. The resulting en- 
tanglement spectrum versus Fermi energy, for both 
choices of partitions, is shown in Figs. |4]3,d. The 
dispersion relation is overlaid on the entanglement 
spectrum as a guide, with the restrictions on A 
schematically emphasized by the solid orange seg- 
ments. The dashed green segments correspond to 
region S, which is traced out. It is clear that the 
left/right mover cut provides the cleaner signature 
(since we are not getting spurious entanglement 
from left /left and right/right coupling) and that 
our interpretation of the delocalized states existing 
where the Fourier components of V n are vanishing 
is correct. 

We have shown that momentum space entangle- 
ment contains valuable information for the identi- 
fication of localized and delocalized states in dis- 
ordered free-fermion models. The next step is to 
generalize the method to higher dimensions to see 
if the entanglement signatures are as clear as in ID. 
It is not obvious whether or not these results will 
easily generalize, but there will be some interest- 
ing features of momentum space entanglement in 
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higher dimensional systems, and especially in topo- 
logical phases with ID edge states and extended 
states buried in a localized bulk state region. This 
method may also provide a nice way to charac- 
terize disordered interacting systems since it only 
requires a ground-state wavefunction and not the 
excited states. Momentum space cuts could be car- 
ried out efficiently for ID interacting systems and 
may yield useful information. 
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